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This is a brief note discussing the energy dependence of superluminal neutrino velocities recently 
claimed by OPERA The analysis is based on the data provided there on this issue, as well as on 

consistency with neutrino data from SN1987a as recorded by the Kamioka detector [3]. It is seen that 
it is quite difficult to reconcile OPERA with SNf987a. The so called Coleman-Glashow dispersion 
relations do not do that well, if applied at all neutrino energies. The so called quantum gravity 
inspired dispersion relations perform far worse. Near OPERA energies both a energy-independent 
velocity, as well as a linear energy dependence with an offset that is comparable in value to the 
observed 5v by OPERA at 28.1 GeV work very well. Our analysis shows that precision arrival 
time data from SN1987a still allow for superluminal behaviour for supernova neutrinos. A smooth 
interpolation is given that reconciles OPERA and SN 1987a quite well. It suggests a fourth power 
energy dependence for ^ of supernova neutrinos. This behaviour is insensitive to whether the 
velocities are energy-independent, or linearly dependent on energy, near OPERA scale of energies. 
Suggestions are made for experimental checks for these relations. 

PACS numbers: 



INTRODUCTION 

The recent announcement by the OPERA collabora- 
tion [ll,!!} that they have measured superluminal neutrino 
velocities has stunned both the scientific community, as 
well as public at large. The announced magnitudes of the 
effects are large and statistically very significant. The full 
implications, should the experiment be vindicated in fu- 
ture, are mind-boggling. Everything in modern physics, 
as we understand it today, depends one way or the other 
on the correctness of the Special Theory of Relativity. In- 
terestingly, speeds greater than c for photons in the same 
energy range as OPERA neutrinos, had been reported by 
Gharibyan [3| before. 

Because of the enormous stakes involved, every ef- 
fort, both experimental and theoretical, must be made 
to make sure that the claims are correct. On the ex- 
perimentalists side, one needs to thoroughly check and 
double check all sources of errors. One needs to repeat 
the experiment, both the original version, as well as vari- 
ants of it like at different energies, different path lengths 
etc.. On the theoretical side efforts have to be made 
to check the consistency of the OPERA data both con- 
ceptually, as well as with earlier precision data. Indeed, 
after the OPERA announcement, a number of interesting 
ideas have been put forward(see [3]-[23|)- Some of these 
ideas even predated the OPERA results [iMl 0, • 

But if the OPERA results are really right, all of physics 
has to be radically revised and there may not exist, at 
present, any sensible way of analysing it. While one may 
come up with some parametrizations, or some pertur- 
bative effective treatments like the one offered by Kost- 
elecky and coworkers jssl-fsij or by Coleman and Glashow 
[iif , the inner selfconsistency of such approaches is never 
guaranteed till another succesful theory replaces Special 
Theory of Relativity. However, theories based on Spon- 



taneous Breaking of Lorentz Invariance 50, 51 [ may be 
better in this regard. That will certainly entail some rad- 
ical revisions of our notions of space and time, to say the 
least. That way we can liken the situation to the days of 
the Bohr atom, where one could parametrize some small 
corrections to observables in the large quantum number 
states, without having the slightest idea of the deep con- 
ceptual changes that the full quantum theory eventually 
brought about. So here too one may have to develop in- 
termediate crutches like the correspondence principle etc 
to make progress. 

To illustrate this point, consider a number of very in- 
teresting papers have appeared which have shown that 
the OPERA result in conjunction with reasonable as- 
sumptions like energy-momentum conservation, can lead 
to contradictory results 52-5^. On the one hand it is 
shown that the kinematics of pion decay is constrained 
to the extent of not permitting neutrinos of the en- 
ergy 'seen' at OPERA, and on the other hand that the 
high energy neutrino looses much of its energy through 
bremsstahlung before it reaches the OPERA detector. In 
these works, some 'reasonable' assumptions like energy- 
momentum conservation, or that charged mesons obey 
normal dispersion relations (based on limits on Cerenkov 
radiation in vacuum), are made. The point is, that in 
the new theory all these notions might themselves un- 
dergo radical revisions. Even if some sort of momentum 
is conserved, it may not be of the form used currently. 
If some sort of relativity principle still survives, that will 
automatically take care of the issue of Cerenkov radia- 
tion in vacuum. In fact such a relativity principle would 
forbid the bremsstrahlung process considered by Cohen 
and Glashow [s^j- According to 57-5l| there are indeed 
Lorentz violating theories that nevertheless are compat- 
ible with a relativity principle. 

In this context, it is worth remembering that the 
OPERA experiment is really probing space-time struc- 
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ture. Though they talk of neutrino energies and such 
other notions, these are not central to an analysis of their 
data. 

Even concepts in which we have great faith like tim- 
ing need to be reassessed. This is essentially a so called 
Time of Flight{TOF) method (is it time to rename it 
as Time of Fright method?), though it is not the tradi- 
tional 'start time - stop time' TOF method, but is instead 
a variant with a start time distribution and a stop time 
distribution. TOF data is mostly analysed using classi- 
cal methods. But the situation here is highly quantum 
mechanical, and it is known that time in quantum me- 
chanics can be a very delicate issue (for various perspec- 
tives on this, see [40|). There are various characterstic 
times like arrival time, passage time etc whose meaning 
in the quantum context can be ambiguous. In the present 
circumstance, neutrinos being very weakly interacting, 
detection time treated quantum mechanically adds fur- 
ther complications. In the end it may well turn out that 
these nuances are not particularly relevant for interpret- 
ing OPERA data, but it is imperative to get a sound 
understanding of these issues given the earth shattering 
nature of the results. 

In a most interesting paper, perhaps the most credible 
interpretation of the OPERA events, Naumov and Nau- 
mov [4lj | have argued that neutrino wavepackets become 
effectively oblate spheroids due to highly relativistic ef- 
fects, and that can lead to a systematic early arrival times 
for neutrinos that are not exactly aligned, at the source, 
with the detector. 

Other subtleties to worry about, regarding the timing 
of neutrinos at OPERA, are conceptual issues tied up 
with the correct way of s ync hronizing clocks 72, IZSl • As 
pointed out by Contaldi [72|, OPERA uses the so called 
one way synchronisation of clocks, which can be par- 
ticularly problematic when synchronizing clocks in non- 
inertial frames such as the earth (which is non-inertial 
both because it is a rotating frame [t^, and also be- 
cause of the earth's gravitational field) . According to this 
analysis, the lack of synchronization due to time dilation 
has the correct sign to explain the OPERA results. But 
building the lack of synchronization to 60 ns from time 
dilation, in the manner described, seems unrealistic, but 
a careful appraisal is certainly called for. 

Before turning to our article, two other proposals, both 
having to do with the beam structure, that could seri- 
ously impact on the interpretation of the OPERA data, 
are worth discussing. The first one, due to Henri [tI ]. 
claims that a fluctuation in neutrino energies of the or- 
der of 10% during the roughly 500 ns. of the leading and 
trailing edges of the PDF can mimic the alleged superlu- 
minal effects claimed by OPERA. While the number of 
protons increase and decrease during these phases, the 
proton energy itself remains at 400 GeV. Thus it is hard 
to see the source of this variation in the neutrino energy. 
Besides, this has to be systematic for the nearly million 



extractions, which is hard to accept. 

The other, potentially more damaging to the OPERA 
interpretation, is the claim by Knobloch [tI] that the 
Beam Current Tracker(BCT) at the CERN end has not 
registered the true PDF, and that even after several sta- 
tistical averaging, 30 and 60 ns. structures are still seen 
in the PDF as measured by the beam current tracker 
(see also 0). If true, this will cast doubts about the 
systematics as claimed by OPERA, and could well nul- 
lify their claims. It is beyond our expertise to assess this 
criticism adequately, and we eagerly await clarification 
by OPERA. 

The plan of this article is as follows: in the next sec- 
tion we recapitulate the most relevant part of the data 
on the supernova 1987a. Here we have only looked at 
the data recorded by the Kamioka detector [4| . We then 
analyse a power law dispersion relation for neutrinos, fix- 
ing the parameters by requiring compatibility with the 
observed differential superluminal neutrino velocity of 
2.48 • 10~5 in units of c at an average energy of 17 GeV 
and the Kamioka data. We then compute the expected 
differentials at the other energies 13.9 GeV, 28.1 GeV 
and 42.9 GeV where OPERA has data. We compare the 
calculated values with the measured values though the 
statistics of these measurements are not as good as at 
17 GeV. We also repeat this exercise by first demanding 
compatibility between OPERA data at 28.1 GeV, and 
then comparing the calculated values with the data at 
13.9 GeV and 42.9 GeV. 

In the section after this, we analyse both a flat, energy- 
independent velocity, as well as a linear energy depen- 
dence taking into account a possible constant shift term. 
In this section we also analyse a linear energy dependence 
with a different shift term which is taken to be valid only 
near OPERA energies. In the subsequent section we dis- 
cuss the difficulties of a treatment that is good both at 
the OPERA energy scales as well as at the SN1987a en- 
ergy scales. We propose, as an illustrative example, a 
form that works at both ends. In the last section we 
make some concluding remarks. 



ESSENTIAL 1987A DATA 

There is a truly enormous amount of literature 
discussing various aspects of the supernova explosion 
SN1987a in the large magellanic cloud in February 1987. 
Fortunately for this analysis we need to concentrate only 
on a few essentials of the detection by Kamioka [3| . We 
briefly summarize them in the next para. 

The distance to the event, Lsn, is 1,68,000 light years. 
A neutrino cluster of 12 events was observed 2 hours be- 
fore the optical sighting ^] of the supernova. The neu- 
trino energies spread over a range 7.5 — 20.0 MeV. In this 
analysis we have excluded the rather low energy event of 
6.3 ±1.7 MeV event arriving 0.686 s. after the burst neu- 
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trino, and three rather high energy events with energies 
35.1 ± 8.0,21.0 ± 4.2,19.8 ± 3.2 MeV, arriving respec- 
tively at 1.541, 1.728, 1.915 s. after the burst event. The 
exclusion of these anomalous events |42l | does not affect 
the reliability of the analysis given here. The first event 
was a Ve event with 20.0 ± 2.9 MeV, identified as a burst 
neutrino. The subsequent neutrinos, the so called ther- 
mal neutrinos comprise of neutrinos and antineutrinos of 
all flavours. The arrival times of the cluster was spread 
around 12 seconds, and this accurately established data 
is crucial to this analysis. Leaving out the anomalous 
events, the highest neutrino energy is around 20 MeV, 
and the lowest around 7.5 MeV. The average neutrino 
energy is around 15 MeV. 

One should not interpret the early arrival of neutrinos 
by two hours compared to photons as evidence for any 
superluminal neutrino velocity as the photons take longer 
to escape the opaque stellar medium. This is similar to 
light taking a million years longer to escape the solar 
interrior. 

The SN1987a neutrinos are of electron flavour, while 
the OPERA neutrinos are predominantly v^. Thus one 
could take the attitude that there is no common ground, 
and hence no real conflict between the two. But os- 
cillations make the situation more complicated, and a 
reasonable attitude to take is that on the average every 
flavour of neutrino would have velocities behaving simi- 
larly in both set ups. There are also indications that the 
OPERA effect has to be by and large flavour indepen- 
dent if serious troubles with oscillation phenomena are 
to be avoided jll]. This provides a justification for the 
analyses to be discussed next. 



Analysis based on photon arrival time 

As mentioned before the first optical sighting of the 
supernova was 2 hrs after the neutrino burst ^]. It is 
possible to use this information to already draw various 
conclusions of interest. 

The analysis proceeds as follows: let iemi^em be the 
emission times, at the supernova, of photons and neu- 
trinos, respectively. Let us for the moment assume that 
photons travel with velocity c. The fact that they actu- 
ally travel slower because of interstellar dispersion, and 
due to a possibly very small mass, will only strengthen 
the conclusions drawn here, not weaken them. Let 
^arr^'^arr bc the rcspcctive arrival times. 

If the neutrinos also travel luminally or subluminally, 
no interesting conclusions can be drawn. However, if the 
SN1987a neutrinos also travelled superluminally, and this 
can not be ruled out, one can use the optical sighting 
data to constrain the superluminality. This will be a 
cruder bound compared to the one to be given in the 
next subsection, but still useful. 

With these assumptions, and an optical sighting 7200 



seconds later, one gets 
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which can be rearranged as 

StsN - 7200 + tl^ - tl^ 
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Since under no circumstances can the photons leave the 
supernova earlier than the neutrinos, one gets 



StsN < 7200 



(3) 



seconds. This is of course being ultra-conservative. One 
may realistically take photons to have left one hour later 
and get 



StsN < 3600 
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168000?/. 



Analysis based on neutrino arrival times 

The analysis now proceeds as follows: let ie^'*, 4*?^ be 
the emission times, at the supernova, of neutrinos of the 
highest and lowest energies, respectively. Let t'^lPj^,t[ 
be the respective arrival times. One now gets 
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(6) 



Unlike the photon arrival case, no simple arguments exist 
for (4"?'' )' ^'^d it would require details of neutrino 

emission mechanisms from supernovae(to get some idea 
of these issues see 43-13]). In most scenarios considered 
in these works, high energy neutrinos precede low energy 
ones. It seems very reasonable to take the quantity [..] 
in eqn.® to be about 10 seconds. This leads to 
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SN 



Sv^^^} < 2.0 



10 
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(7) 



It is also important to consider the 1MB data [6| in a 
more careful analysis. The energies of the detected neu- 
trinos are higher in 1MB. If very accurate timing was 
available for all detectors, that also could have been used 
to constrain the energy dependence of Svsn- But unfor- 
tunately, the Kamioka timing had an uncertainty of ±1 
min; the 1MB had by far the most accurate timing record, 
with an error of ±50 milliseconds. A determination of the 
so called offset time from fitting data to models of neu- 
trino emission seems to indicate that the first 1MB neu- 
trino arrived 0.1 s. before the first Kamioka neutrino, in 
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the so called two component model of neutrino emission 
(47} . But given the complexity of the issues, too much 
significance should not be attached to this. An intriguing 
suggestion has been made by Fargion [s^ to search past 
records of 1MB for evidence supporting superluminal be- 
haviour of neutrinos. 

If a more detailed analysis of arrival times could be 
made on the basis of a soundly established model of neu- 
trino emission from the supernova, one could put even 
more stringent limits than the one in eqn.Q, but the 
situation here is very complex (see for example [431-147']). 

It should be emphasized that both eqn.(l4|) and eqn.(l7]) 
only reflect the possible upper limits to Svsn- They do 
not imply that supernova neutrinos are superluminal, but 
that the available data are indeed consistent with them 
being so. They may in the end turn out to be subluminal, 
or just luminal too. We have left out any discussion of the 
effect of neutrino masses. The current thinking on these 
masses is that they can at most be of a few eV with one 
of them even being massless. For 20 MeV neutrinos, the 
correction to — is ~ ^ ■ 10~^*, if the mass is taken to 
be 2eV. 



GENERAL POWER LAW 

Now we consider a possible power law dependence for 
6v. This analysis should be taken as an exercise in phe- 
nomenology with no pretenses to any deeper theoretical 
considerations. For integer powers, these belong to the 
class of so called Distorted Special Relativistic theories 
introduced in 57 -HI]. OPERA data has been analysed 
in the context of these models in [i^. An interesting 
point stressed in [4§| is that these DSR type theories still 



allow for a relativity principle. 

The simplest candidate for a power law behaviour is 



— — ai E 

c 



(8) 



where ai is a constant with dimensions of i? We shall 
require equ.® to be valid at all energies. Therefore 



Estimating qi, m from photon arrival time data 

We shall estimate the parameters ai,m first using 
the photon arrival time data. Taking average energy of 
SN1987a at 15 Mev and the OPERA energy at 17 GeV, 
where the statistics are the best, the energy ratio comes 
out to be 8.8 • 10^^. Taking the estimate of eqn.(|4]) for 
SvsN one gets, on using eqn.®. 



TO = 1.5 ai^(20TeV) 



-1.5 



(10) 



Using eqn. pH]) . one can estimate S(3 = — , in units of 
10"^ at the energies of 13.9, 28.1, 42.9 GeV : 

(5^(13.9) = 1.8 (5^(28.1) = 5.2 5^(42.9) = 10.0 

(11) 

respectively. Thus the estimates based on the powerlaw 
fit of eqn.® are off the experimentally measured cen- 
tral values by .5, 9, 10 standard deviations(statistical) 
respectively. Thus this power law fit works rather poorly. 

This can be understood from the fact that the val- 
ues at 17GeV and 28.1GeV are very close. In fact, they 
are consistent with no energy dependence at all. It is 
therefore instructive to make an analysis with the data 
at 13.9, 28.1, 42.9 GeV only. Of course the statistical 
errors are much higher, and this may look like deliber- 
ate data selection, but its worth going through it given 
the scantiness of data at different energies. An analysis 
identical to the one above, now gives 



m = 1.39 ai = {50TeV) 



-1.5 



(12) 



Using eqn. p^ . one can estimate (5/3 = —, in units of 
10-^ at the energies of 13.9, 42.9 GeV : 



(5/3(13.9) = .93 (5/3(42.9) ^ 4.46 



(13) 



respectively. Now the estimates based on the powerlaw 
fit of eqn.® are off the experimentally measured cen- 
tral values by 1.63, 2.3 standard deviations (statistical) 
respectively. As a last bit of comparison, we estimate 5/3 
at llGeV that eqn.iO would give: (5/3(17) = 1.39- 10-^ 
differing by the measured value by 3.9 standard devia- 
tions. 



Svop Eqp 



(9) 



Estimating ai, m from neutrino arrival time data 



The data from OPERA for the advancement 
times are: 60.7, 60.3, 53.1, 67.1 nanoseconds at 
17.0, 28.1, 13.9, 42.9 GeV respectively. In terms of 
^ these correspond to (2.48, 2.46, 2.17, 2.74) • 10"^ 
respectively. The statistical errors are the least 
at 17 GeV, getting progressively higher at the 
other energies. The quoted statistical errors are : 
(±0.28, ±0.54, ±0.77, ±0.75) • lO'^ respectively, and 
the systematic errors arc ± 0.30 • 10~^ [l|. 



In a similar manner we can estimate the parameters 
ai, m by using the estimates based on neutrino arrival 
times as given in equ.©. As input, we need Svsn at the 
highest energy, 20 MeV. Getting it from eqn.dT]) requires 
an estimate for Svsn at the lowest energy, taken here 
to be 7.5 MeV, which can be obtained as a fraction of 
the former if we already know the exponent m. We shall 
circumvent this in a practical manner by simply taking 
the rhs of eqn.© to be the estimate for (5/3 at 20 MeV. 
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As before, using this estimate for S[3gj^ in eqn.® gives 

m = 2.42 ai = {lATeVy^-*^ (14) 

Using eqn. (fT4|) . one can estimate 3/3 — in units of 
10"^ at the energies of 13.9, 28.1, 42.9 GeV : 

<5/3(13.9) = 1.52 (5/3(28.1) = 8.4 ,5/3(42.9) = 23.3 

(15) 

respectively. 

Thus the neutrino arrival times give a much steeper 
power law than the one given by photon arrival times. 
It is therefore not surprising that the estimated results 
are much worse compared to those given by eqn. (jll[) . As 
before, we give the results that one would get if one had 
only used the data at 13.9, 28.1, 42.9 GeV: 



m = 2.11 ai = {AATeV) 



-2.11 



(16) 



Using eqn. pB]) . one can estimate 6(3 — —, in units of 
10-^ at the energies of 13.9, 42.9 GeV : 



(5/3(13.9) = .56 (5/3(42.9) = 6.00 



(17) 



Quadratic corrections of the type m — 2 have been sug- 
gested countless number of times in the so called Quan- 
tum Gravity inspired models of dispersion relations. If 
OPERA data even at high energies is taken seriously, 
these are seen to perform rather poorly. Furthermore, 
the intrinsic energy scales in these models is the Planck 
Mass, but in these fits that scale is turning out to be 
around TeV. But then, there are also models where the 
'Planck Mass' can be much lower. 

The lesson one learns from the analysis so far is that 
while fitting SN1987a requires steeper energy depen- 
dence, the OPERA data on energy dependence prefers 
flatter dependence. Striking a balance between these op- 
posing trends is the key to reconciling the OPERA and 
SN1987a data. 



FLATTER ENERGY DEPENDENCES. 
Energy Independent Superluminality. 

If one takes both the data of 60.7 nanoseconds at 17 
GeV, as well as the 60.3 nanoseconds at 28.1 GeV, one 
has to more or less conclude that the observed superlu- 
minality is energy independent: 



6(3 



OP 



2.48 • IQ- 



(18) 



This is what the dispersion relations from the so 
called Standard Model Extensions 35l-l37l| as well as the 
Coleman-Glashow analysis [i^ would give. These pro- 
posals were made in the context of the so called High 
Energy Violations of Lorentz invariance. But for neutri- 
nos, these analyses would seem to imply their validity as 



long as neutrino energies are much higher than neutrino 
masses. Hence they should be applicable to the SN1987a 
neutrinos too. 

A flat velocity surplus of eqn.(|18|) would imply that 
the supernova neutrinos preceded the photons by about 
4 years! While neutrinos hitting the earth 4 years before 
the photons can not be ruled out the fact that the 
neutrinos observed by Kamioka almost fully accounted 
for the energy output in neutrinos that is expected from 
the standard supernova models, would cause some prob- 
lems of energetics. But it is fair to say that one needs 
a complete analysis of neutrino oscillations with the su- 
perluminality of eqn. (|18|) built in before one can make 
confident statements. There are already such analyses of 
neutrino oscillations available in the literature 37, 6^ . 
But this line of thinking requires the above equation to 
be valid at all energies. But we shall argue that it is es- 
tablished, if at all, only at OPERA energies, and we shall 
propose interpolations that can comfortably sit with the 
SN1987a data presented earlier. 



Linear Energy Dependences. 

It is interesting that the central values of 6t quoted 
by OPERA at 13.9 and 42.9 GeV, namely, 53.1 ns and 
67.1 ns, when combined with 60.3 ns at 28.1 GeV, sug- 
gest a rather good linear fit to the data. But from our 
analysis till now, a fiatter energy dependence at OPERA 
energy scales does not extrapolate well to the supernova 
scales. In fact the linear dependence suggested by the 
above mentioned data is roughly 1 ns increase for every 
2 GeV. Extrapolating 53.1 ns at 13.9 GeV to supernova 
energy scales (nearly zero on OPERA energy scales) , one 
still gets 46 ns which translates to a ^ of ~ 2 • 10~^. 
This is unacceptably large, leading to neutrinos arriving 
nearly 3.4 years before photons(see, however, f68!|)! 

A simple way of circumventing this is to consider a 
linear dependence with a suitable offset 



OP 



a2E + 6(3q 



(19) 



where 6(3q is so adjusted as to cancel most of the large 
value mentioned before. In eqn. p^ . ai is roughly 4.0 • 
10~^ per 2 GeV. Even without working out (5/3q, we can 
see that this too ends up in trouble. This is because 



c ohinh 



r alow 



ai ■ 12.5MeV ~ 2.5 • 10 
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(20) 



is still very large, predicting a differential time of arrival 
of nearly two hours for these neutrinos, whereas Kamioka 
data constrains them to be at most 12 seconds. 

Eqn. (fT9)) . with just the offset term, is what the 
Coleman-Glashow like Lorentz violating dispersion rela- 
tion 35, H, 13] would give. Thus unless one interprets 
the Coleman-Glashow form to be valid only for energies 
well above a GeV scale, we see that it does not work for 
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neutrinos. Egn.dTOl 



is an example of the so called DSR 
Among the many 



type theories discussed by |57H5i 
Lorentz violating theories, these may have the advantage 
that they may still accommodate a relativity principle. 

As it appears to be very difhcult to reconcile both 
OPERA and SN1987a data with a single scale energy 
dependence, let us abandon that approach and first seek 
an energy dependence valid close to the OPERA ener- 
gies (namely 10 to 40 GeV) only. Then, a very nice 
fit to the available data is provided by eqn.(fT9|) with 
(5 ■ lO^TeV)-^ and ^ = 1.91 ■ 10-^ that is 



a2 



6/3op = 2-10 
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E 
GeV 



1.91 • 10" 



(21) 



A comparison of OPERA data with various forms of 
dispersion relations has also been done in [4^. Their 
comparison with SN1987a is rather preliminary, and not 
as detailed as the one given here. They even hint at some 
sort of a metatheorem{\) claiming that it would not be 
possible to achieve consistency with both OPERA and 
SN1987a data. We do not agree with that. Many pa- 
pers have appeared that have analyzed the compatibility 
issues between OPERA and SN1987a ^-M]- 



INTERPOLATING OPERA AND SN1987A 

Thus, if we want to give an energy dependence formula, 
with a single intrinsic scale, that fits both OPERA and 
SN1987a data, it is almost impossible. If , however, one 
allows another intrinsic scale, many possibilities open up. 
We do not wish to explore all possible such forms, as such 
an exercise is somewhat pointless at the present juncture. 
We give one of the simplest forms: 



^GeV 



^GeV 



Ml 



(22) 



where S/3qp is from either ean. (IT5|) or ean. (pT|) . and Mj^ 
in GeV is a scale which is much larger than supernova 
energies, but much smaller than OPERA energies. The 
parameters Ml, n are constrained by SN1987a data ac- 
cording to 
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SN 



p2 

2.0- 10-5 (-2^)" 



(23) 



We have used ~ here because the prefactor in eqn. (P5)) 
is 2.48 if we had used eqn.([TS|. and 1.91 if we had used 
egn. d^IT) . This hardly makes much difference to the esti- 
mates of Ml, n. We will see in the next subsections that 
estimates of eqn.Q and equ.© quite severely constrain 
the parameters in eqn. 



Photon arrival constraints 

Let us begin with the estimate as given in eqn.Q and 
take it to be for 15 MeV neutrinos. It is straightforward 



to work out that for n = 1,2 the scale Ml works out 
to 2.5, 0.2 GeV respectively. Both these scales satisfy 
our criterion that the scale should be much smaller than 
OPERA energies, and much larger than supernova neu- 
trino energies. At this stage there is no other criterion to 
choose one over another. 



Neutrino arrival constraints 

If instead of eqn.(|3]) we use eqn. ([7]), we get that for 
?i = 1, 2 the scale AIl is 50, 1 GeV respectively. Clearly, 
Ml = 50 GeV is unacceptable, while Ml — 1 GeV is 
certainly acceptable. In any case, one should be using the 
more accurate and cleaner estimates provided by eqn.©, 
and rule out n 1 in eqn. 



DISCUSSION AND CONCLUSIONS 

Our aim was to discuss possible energy dependence 
of the superluminal neutrino velocities reported by the 
OPERA collaboration. Our motivations are three fold: 
i. an understanding of energy dependence will provide 
an important key to unraveling the mysteries, ii. without 
any energy dependence it is impossible to reconcile the 
findings of OPERA with the accurate timing data from 
SN1987a, and finally, iii. if OPERA data at both 17 and 
28.1 GeV are taken seriously, then, near OPERA energies 
the velocities are energy-independent. But the mean val- 
ues they quote for 13.9, 28.1 and 42.9 GeV hint at a clear 
energy dependence. Though the statistical errors for the 
measurements at 13.9 and 42.9 GeV are much poorer, 
as the sample sizes were drastically reduced in number 
due to the reliance on only CC events, we have taken 
the central values, which are significant at three stan- 
dard deviations, to parametrize the energy dependence 
on the one hand, and to reconcile them with SN1987a on 
the other. 

For the SN1987a data, we have used both the timing 
of the first optical sighting, as well as the time spread in 
the cluster of 12 events observed at the Kamioka detec- 
tor. We have first attempted a power law fit to both 
the OPERA data at 17 GeV as well as the SN1987a 
estimates. Fitting the light arrival time gives an in- 
trinsic mass scale of 20TeV and an energy exponent 
of 1.50. The calculated values at 13.9, 28.1 and 42.9 
GeV from this fit yield values that deviate from the 
quoted mean values by .5, 5.2 and 10 standard devia- 
tions (statistical). On the other hand, fitting photon ar- 
rival times to OPERA data at 28.1 GeV yields a scale of 
50 TeV, and an exponent of 1.39. 

A fit based on the 12 seconds clustering of the neutrino 
events at Kamioka iand 17 GeV data from OPERA gives, 
on the other hand, an intrinsic mass scale of 1.4 TeV and 
an exponent of 2.42, which could be interpreted as really 
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2.0 given the tentativeness of the estimates. Extrapo- 
lating these to 13.9, 28.1 and 42.9 GeV, one finds the 
calculated values of 5f3 to be 1.52, 8.4 and 23.3. Such 
dependences with m = 2 are suggested by many quan- 
tum gravity inspired models, but the scales found here 
are very different. 

Motivated by the near linear increase of 5t reported 
by OPERA for the observations at 13.9, 28.1 and 42.9, 
we next attempted a linear fit but with an offset so that 
the linear behaviour extrapolated to supernova energies 
would reconcile with the supernova arrival limits. If the 
slope is determined by the OPERA data, it is found that 
irrespective of the offset, one gets a differential arrival 
time for 7.5 MeV and 20.0 MeV supernova neutrinos that 
is very clearly in conflict with the arrival times. 

Interpreting all this as the inability of any single-scale 
energy dependence to account for both OPERA and 
SN1987a data, we decided to check how well a linear 
fit will work only when restricted to OPERA energies. 
A mass scale of 5 • 10^ TeV and an intercept (offset) of 
2.0 • 10^^ was obtained. This fit works very well. This 
can be tested by measuring 5t at other, but comparable, 
energies as OPERA. 

Finally, an interpolation was sought between this lin- 
ear behaviour valid near OPERA energies, and the up- 
per limits on neutrino velocities obtained from supernova 
data. The important point to notice is that these super- 
nova arrival time limits are still consistent with superlu- 
minal velocities for supernova neutrinos. This is an ex- 
citing possibility that needs to probed further. There are 
clearly many such interpolations possible. One of them 
is given in ean. (l^51) . The additional scale Ml and the 
exponent n are again constrained by arrival time data 
and the following. 



5/3 



-P2 

^GeV 



1.0' 



(24) 



seems most probable. Thus, this analysis anticipates an 
E"^- dependence for 5v near supernova energies. Even at 
1 GeV, this relation predicts a — which is only 
that reported by OPERA at 28.fGeV! 



ith of 
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Further Work 

The OPERA collaboration should increase the statis- 
tics of the measurements done with the CC events. De- 
termining reliably the energy dependence of the effect, 
if true, is one of the most important things to be done. 
Efforts should also be made to measure the effect at dif- 
ferent energies. An accurate determination of the energy 
distribution of the neutrino flux is also important. 

It is important also to perform this experiment at dif- 
ferent baseline lengths because it is crucial to separate 
effects of production and detection from those of propa- 
gation(the effects discussed by 74| are an example of the 



kind that have no dependence on the baseline distance). 
The claim of 18cm accuracy in the baseline length is cen- 
tral to the establishment of the superluminality. It is 
therefore critical to flnd several independent checks on 
what could easily become the Achilles heel for the exper- 
iment. 

As already emphasized by Coleman and Glashow [il] , 
the superluminal velocity effects and the energy depen- 
dence of them can have novel manifestations in oscilla- 
tion phenomena. It is therefore of utmost importance 
to remodel oscillation calculations by including the en- 
ergy dependences considered here, and then subject them 
to experimental tests. In principle, the various param- 
eters that entered the interpolation formula in eqn. (j23p 
could be flavour dependent. Even the TOP analysis itself 
should be remodelled taking into account possible flavour 
dependences, as well as effects of oscillations. 

The supernova neutrino data, from all the detectors 
that observed them, also needs to be thoroughly reexam- 
ined, particularly from the points of view of energetics 
and timing. 

The most challenging task, should OPERA results be 
found to stand further scrutiny, will be in arriving at a 
reappraisal of our notions of space and time. The entire 
edifice of theoretical physics stands on these pillars. 
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